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A simultaneous description of non-strange nuclei and hypernuclei is provided by a single mass 
formula inspired by the spin-flavour SU(6) symmetry breaking. This formula is used to estimate the 
hyperon binding energies of Lambda, double Lambda, Sigma, Cascade and Theta hypernuclei. The 
results are found to be in good agreement with the available experimental data on 'bound' nuclei 
and relativistic as well as quark mean field calculations. This mass formula is useful to estimate 
binding energies over a wide range of masses including the light mass nuclei. It is not applicable for 
(-H ■ repulsive potential. 
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The hypernuclear physics is of great importance in many branches of physics. Of particular interest is the un- 
derstanding of strange particles in baryonic matter, since many questions in heavy-ion physics, particle physics and 
astrophysics are related to the effect of strangeness (S) in nuclear matter. Moreover, the contribution of the hyper- 
ons strongly influences the mass of neutron stars as well. In the past decade considerable amount of spectroscopic 
informations were accumulated experimentally on the A (S=-l)hypernuclei. The A separation energies were deter- 
mined for the ground states of about 40 A hypernuclei including several double- A hypernuclei [1], [2]. Doubly strange 
hypernuclei also arise in a form of S~ (S=-2) hypernuclei and were studied both theoretically [3] and experimentally 
[1], [4], [5] in a limited number of nuclei. Sigma hyperons exhibit interesting property in its interaction with nuclei. 
Several studies suggest that due to strongly repulsive E-nucleus potential Sigmas are unbound in nuclei, except for 
the very special case of nuclei with mass number A=4 [6]. Existance of a T= 1/2, S= bound state in %He was first 
predicted by Harada et al. [7] and has been found experimentally [8], [9], [10], [11]. The 1/A dependence of the Lane 
term present in the isospin-dependent S-nucleus potential [7] reduces the possibility of finding bound E-hypernuclei 
with large A. 

In this scenario, an exotic hyperon + (S=+l, mass ~1530 MeV, width <15 MeV ) with exotic pentaquark structure 
was predicted in 1997 [12]. Announcement of its discovery at Spring-8, Japan [13] sparked an avalanche of activities 
in the field of hyperons and hypernuclei [14]. Now a question has arisen whether the + hyperon really exists [15] 
and if so, whether it will be bound in a nucleus [16], [17]. Calculations in a relativistic mean-field formalism (RMF) 
suggest that as there is an attractive + -nucleus interaction, the + particle can be bound in nuclei and, the + 
hypernuclei would be bound more strongly than A hypernuclei [18]. Recent calculations in quark mean-field (QMF) 
model also support existence of bound + hypernuclei and predict that in comparison to A hypernuclei more bound 
states are there in + hypernuclei [19]. While a search for bound Theta hypernuclei is on, for a large number of 
hypernuclei, including double Lambdas, Cascade and Sigmas, more experimental data are needed. One also needs to 
have an apriori estimation of their possible binding energies in a wide mass range for planning of the experiment and 
to locate the peaks in the experimental missing mass spectra. 

Although the A, S, 3 and O hyperons are all baryons, no single mass formula exists which can predict the binding 
energies of all of them on the same footing as the non-strange nuclei. In this work we present a generalized mass 
formula for both the strange and non-strange nuclei which does not disregard the normal nuclear matter properties. 
It is a straight forward equation and binding energy of non-strange as well as bound hypernuclei can be estimated 
by using this single equation. Its predictions compare well with the available experimental data. However it is not 
applicable for repulsive potential. 

Earlier, Dover et al [3] prescribed two separate mass formulae for A and S hypernuclei by introducing several volume 
and symmetry terms in Bethe-Weizsacker mass formula (BW). The BW formulae given in Ref. [3] were developed 
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for a much broader context (multi-Lambda systems or strange hadronic matter) with parameters inspired on the 
strengths of the YN and YY interactions. In this respect, the single-Lambda BW equation given in [3] is a limit of 
a more general formula. The BW equation proposed in Ref. [20] is inspired by the spin-flavour SU(6) symmetry and 
the pairing term is replaced by the expectation value of the space-exchange or, Majorana operator. The Majorana 
operator in the ground state configurations of single- and double-Lambda hypernuclei were expressed in terms of N, 
Z, and Lambda number. A strangeness dependent symmetry breaking term was also incorporated. This prescription 
gave reasonable description of the experimental data on singlc-and double-Lambda hypernuclei separation energies. 
As none of these formulations had explicit hyperon mass consideration, they can not be used for binding energy 
calculation of other hypernuclei. Both the formulae have some other limitations which will be discussed later. 

The Wigner's SU(4) symmetry arises as a result of the combined invariance in spin (I) and isospin (T). In order 
to incorporate the strangeness degree of isospin, SU~t(2) is replaced by SUp(3) and the combined spin(I)-flavour(F) 
invariance gives rise to the SU(6) classification of Gursey and Radicati [21]. The SUf(3) symmetry breaks by explicit 
consideration of a mass dependent term in a mass formula. The SU(6) symmetry breaking is related to different 
strengths of the nucleon-nucleus and hyperon-nuclear interactions and has important consequences. For example, 
although small, the X — A mass difference figures prominently in the smallness of the A-nuclcar spin-orbit interaction 
[22]. 

In this work the non-strange normal nuclei and strange hypernuclei are treated on the same footing with due 
consideration to SU(6) symmetry breaking. The generalization of the mass formula is pursued starting from the 
modificd-Bcthc-Wcizsacker mass formula (BWM) preserving the normal nuclear matter properties. The BWM is 
basically the Bethe-Weizsacker mass formula extended for light nuclei [23], [24], [25], [26] which can explain the gross 
properties of binding energy versus nucleon number curves of all non- strange normal nuclei from Z=3 to Z=83. In 
BWM, the binding energy of a nucleus of mass number A and total charge Z is defined as 

B(A, Z) = a v A - a s A 2 ' 3 - a c Z(Z - I)/ A 1 ' 3 - a sym (N - Z) 2 /[(l + e- A / k )A] + 5 new , (1) 
where for normal nuclei N and Z are the number of neutrons and protons respectively and 

a v = 15.777 MeV, a s = 18.34 MeV, a c = 0.71 MeV, a sym = 23.21 MeV and k = 17, (2) 
and the pairing term, 

Snew = (1 - e- A/c )S, where c = 30, (3) 

and 

S = YlA^ 1 !" 1 for even neutron — even proton number, = — Y2A~ X I 2 for odd neutron — odd proton number, 
= when total neutron plus proton number is odd. 

(4) 

Hypernuclei are found to be more bound than normal nuclei. A systematic search of experimental data of hyperon 
separation energy (Sy) for A, AA, E° and S~ hypernuclei leads to a generalised mass formula for hyper and non- 
strange nuclei which will be, henceforth, called the BWMH. The experimental Sy of the A-hypernuclei (for which 
experimental data are available over a wide mass range) is found to follow a relation Sy cx A~ 2 I 3 [27]. The available 
data for Cascade hypernuclei also follow a similar Sy cx A~ 2 I 3 trend but with different slope. Unlike Levai et al. 
[20], the SU(6) symmetry breaking term represented by the strangeness is taken here to be inversely proportional to 
A 2 / 3 . Explicit inclusion of the pairing term partly accounts for the Majorana term while still preserving the nuclear 
saturation properties. An additional mass dependent term breaks the SUf(3) symmetry. 

In BWMH the hypernucleus is considered as a core of normal nucleus plus the hyperon(s) and the binding energy 
is defined as 

B(A, Z) = a v A - a s A 2 ' 3 - a c Z{Z - I)/ A 1 ' 3 - a sym (N - Z c ) 2 /[{I + e- A ' k )A] + S new 
+n Y [c .(m Y ) - ci - c 2 | S \ /A 2 / 3 }, 

(5) 

where n Y — number of hyperons in a nucleus, m Y — mass of the hyperon in MeV, S = strangeness of the hyperon 
and mass number A = N + Z c + n Y is equal to the total number of baryons. N and Z c are the number of neutrons 
and protons respectively while the Z in eqn.(5) is given by 

Z = Z C + n Y q (6) 
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where q is the charge number (with proper sign) of hyperon(s) constituting the hypernucleus. For non-strange (S=0) 
normal nuclei, Z c = Z as ny =0. The choice of S new value depends on the number of neutrons and protons in both 
normal and hypernuclei. For example, in case of \Li the neutron number N=5(odd), proton number Z c =3(odd), 
and ny = l. Therefore, 8 = — as the (N, Z c ) combination is odd-odd, although the total baryon number A= 

A = N + Z c + ny =9(odd). Whereas, for non-strange normal 9 Li nucleus 5 = for A=9(odd). 

In eqn.(5), the constants Co,c\ and C2 have been fixed from an empirical fit to the experimental values of Sy for 
thirty-five A, three A — A and six S~ hypernuclei (Table 1). 



TABLE I. 
of hyperons 


Experimental data on hyperon separation energies (Sy) of hypernuclei 
(ny) and NS= ny-S where S is the strangeness. 


with experimental errors 


(ASy), number 


Hypernuclei 


Sy 
MeV 


ASy n Y 
MeV 


NS 


Ref. 


\H 


2.04 


0.04 1 


-1 


[1] 


\He 


2.39 


0.03 1 


-1 


[1] 


\He 


3.12 


0.02 1 


-1 


[1] 


%He 


4.18 


0.10 1 


-1 


[1] 


\He 


5.23 


0.00 1 


-1 


[28] 


iHe 


7.16 


0.70 1 


-1 


[1] 


iLi 


4.50 


0.00 1 


-1 


[1] 


7 A Li 


5.58 


0.03 1 


-1 


[1] 


iLi 


6.80 


0.03 1 


-1 


[1] 


9 A Li 


8.50 


0.12 1 


-1 


[1] 


r A Be 


5.16 


0.08 1 


-1 


[1] 


iBe 


6.84 


0.05 1 


-1 


[1] 


iBe 


6.71 


0.04 1 


-1 


[1] 


i°Be 


9.11 


0.22 1 


-1 


[1] 


iB 


8.29 


0.18 1 


-1 


[1] 


i°B 


8.89 


0.12 1 


-1 


[1] 


i'B 


10.24 


0.05 1 


-1 


[1] 


i 2 B 


11.37 


0.06 1 


-1 


[1] 


l 2 c 


10.76 


0.19 1 


-1 


[1] 


A ^ 


11.69 


0.12 1 


-1 


[1] 


lA r < 
A 


12.17 


0.33 1 


-1 


[1] 




12.17 


0.00 1 


-1 


[1] 


i 5 N 


13.59 


0.15 1 


-1 


[1] 
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fo 

A 


12.50 


0.35 




-1 


[291 


A ^ 


13.59, 


0.0 




-1 


[301 


a Si 


16.00 


0.29 


1 


-1 


[291 

L J 


fs 


17.50 


0.50 


1 


-1 


[11 

L J 


A 


17.96 


0.00 


1 


-1 


[301 

L J 


A °Ca 

A w 


18.70 


1.10 


1 


-1 


[291 

L 


A CO 


19.24 


0.00 


1 


-1 


[301 


A * 


19.90 


1.00 




-1 


[311 

L J 


5 A 6 fe 

A 


21.00 


1.50 




-1 


[311 

J 


fY 


22.10 


1.60 


1 


-1 


[311 

L w J 


\ 39 La 

A j-"-" 


23.8 


1.00 


l 


-1 


[321 

L" J 


| 08 Pfe 


26.5 


0.5 




-1 


[321 

L w J 




7.25 


0.19 


2 


-2 


[2] 


'"Be 

A A ^ 


17.7 


0.4 


2 


-2 


[11 

L J 


13 D 

A A 


27.5 


0.7 


2 


-2 


[33,34] 




5.90 


1.2 




-2 


[11 




9.2 


2.2 




-2 


[11 

L J 


B-C 


18.1 


3.2 




-2 


[il 

L J 


15 ^ 


16.0 


4.7 




-2 


[i] 


B-C 


16.0 


5.5 




-2 


[i] 


a 8 _AZ 


23.2 


6.8 




-2 


[i] 



The convention used in Table 1 (as well as in the text) is y q Z where Z is the net charge, Y is the hyperon type 
with charge q and, A is the total number of baryons [1,7]. There exists another convention in the literature in which 
one uses y q Z c for denoting the hypernuclei, where Z c is the number of protons. Therefore, according to this other 
convention, the hypernucleus ^+He of the present manuscript will be ^+H. Also, q+Li would be Q+He, and Q+Be 
would be Q + Li. 
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Here we calculate the root mean square deviation i.e., r.m.s. (<r) for the hyperon separation energies where, 
a 2 = (l/iV)S [(SV)t?i. — (Sy)ex.] 2 ■ Due to small number of data, simultaneous 3 parameter search does not yield 
any meaningful result. After many trial searches we fixed the value of cq = 0.0335 and made two parameter search 
using 35 A hypernuclei data which yielded ci=27.59 (45) and ci= 48.6 (21) with r.m.s. =1.416 MeV. For better results 
from actual plotting of the graphs, ci=26.7 was chosen and, keeping it fixed, a two parameter search with 35 A 
hypernuclei data lead to co=0.0327 (4) and C2=48.6 (21) with r.m.s. deviation 1.416 MeV. Finally we fixed ci=26.7 
and C2=48.7 and performed a one parameter search with 35 A hypernuclei data. This leads to co=0.0327 (2) with 
r.m.s. =1.375 MeV. We fixed co=0.0335 which yielded best results for the two parameter fit. Finally, from the overall 
best fit to all the data (r.m.s ~ 1.4 MeV), we choose co=0.0335; ci=26.7 and C2=48.7 in eqn.(5). 

B(A, Z) = a v A - a s A 2 ' 3 - a c Z{Z - I)/ A 1 ' 3 - a sym (N - Z c f /[(l + e- A ' k )A] + S new 
+n y [0.0335my - 26.7 - 48.7 | S \ /A 2/ % 

(7) 

The hyperon separation energy Sy defined as 

Sy = B(A, Z)h yP er — B(A — fly, Z c ) core , (8) 

is the difference between the binding energy of a hypernucleus and the binding energy of its non-strange core nucleus. 

It is interesting to note that this single equation yields the values of Sy in reasonable agreement with the available 
experimental data of all known bound hypernuclei. Fig.[l] shows plots of Sy versus A for A and AA hypernuclei 
which are in good agreement with the experimental data [1], [2]. The recent discovery of the \°Li nucleus points to a 
value of Sy 10 - 12MeV [35] where the BWMH predicts S A [ A °Li] = 11 AM eV. The BWMH prediction of Sy for 
the AA He also agrees well [Fig. 1(b)] with the recent experimental value 7.3 ± 0.2MeV of Takahashi ct al. [2]. 

Harada et al. investigated the structure of \Be by the coupled-channel calculation between the 3 H + S + and the 
3 He + E° channels and predicted the E + binding energy in fjife to be 3.7 to 4.6MeV [7]. Available experimental 
binding energy values are 4.4 ± 0.3 ± lMeV [8], 2.8 ± 0.7MeV [9], 4 ± lMeV [10]. The BWMH predicts binding 
energy of S° (my = 1192.55MeV) and £+ (my = 1189.37MeF) in f, i2e (=S° + \He) as 2.69MeV and 
(=S + + \H) as l.QMeV respectively. Search for heavier S hypernuclei has been carried out by several authors [6], 
[8], [9], [10], [11], [36] without success. P. K. Saha et al. studied inclusive (ir~ , K + ) spectra on C, Si, Ni, In and 
Bi targets and concluded that a E-nucleus potential is strongly repulsive in such relatively heavy nuclei [6]. It has 
been suggested that the binding of the Sigma in the light hypcrnuclear species like %He is mainly due to a strong 
isovector term of the Sigma-nucleus potential which emphasizes the differences of the Sigma-N interaction in the 
isospin T=l/2 and the T=3/2 channels. This term is not present in the BWMH formula. Nevertheless, without 
changing any parameter, this single general equation (BWMH) reproduces the binding energy of the 'bound' |-ffe 
hypernucleus. As mentioned before, the theoretical formulation suggests that the binding energy of the light Sigma 
hypernuclear species is provided by the isospin dependent Lane term (associated to the isospin asymmetry of the 
core) not included in the BWMH formula. Thus the binding energy of BWMH obtained for light hypernuclei comes 
from the isospin-averaged (i.e. isospin independent) term which provides attraction, while recent analysis on heavier 
Sigma hypernuclei (dominated by this isospin independent term of the Sigma-nucleus potential) suggest repulsion. For 
S° + 2 H, £+ + 2 H, and XT + 2 H hypernuclei the separation energies predicted by BWMH are -3.53 MeV, -4.62 MeV 
and -3.37 MeV respectively indicating that these light Sigma hypernuclei would be unbound, even if the potential is 
attractive. In fact, so far no bound state of these hypernuclei could be found in the experiment. Indeed, further data 
on Sigma hypernuclei are necessary to determine more conclusively whether the Sigma feels attraction or repulsion. 

Recent experimental data for S~ + 12 C and S~ + n B are not conclusive and predict values in the range of 3 to 
10 MeV depending on the S~ well depth [4] while, Khaustov et al. suggests potential depth to be less than 20 MeV 
[5] . In absence of any new conclusive data, we used the available data tabulated in Ref . [1] . BWMH estimations for 
the S - separation energies compare well with the experimental values (Fig. 2a). As no experimental data exists so 
far on the bound Theta hypernuclei, the 9 + separation energies are compared with the recent theoretical predictions 
and the Sy of 9+ are found to be close to the quark mean field (QMF) calculations [19]. 

Analysis of each term of eqn.(8) reveals that the 5 term has very little contribution (positive) in Sy. The a sym term 
difference also contributes positive but, quite small except at very high | (N — Z c ) |. The a s term difference (asd) 
is always negative and being more so at small A (Fig. 3). For q=0, the Coulomb term difference (acd) contributes 
positive but very small for all A. For q=-l it is positive and increases rapidly with A while for q=+l, the acd is 
negative and decreases rapidly with A. In the 6 + plot [Fig. 2(b)] a maximum at low A and decrease in Sy values 
at large A arise mainly due to increasing contribution from the "strangeness term" (i.e., sterml= —48.7 | S \ /A 2 / 3 ), 
surface term (asd) and decreasing (negative) contribution from the Coulomb term (fig. 3c). For q = and —1 such 
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maxima in Lambda plot [Fig. 1(a), 1(b)] and S~ plot [Fig. 2(a)] are absent as asd, acd and sterml terms make 
increasing contributions with A (fig. 3a, 3b). 

Interestingly, the large masses of hyperons make some strongly bound hypcrnuclei whose non-strange core nuclei 
might be unbound. For example, BWMH shows that the prediction for the hyperon separation energy Sy for the 
Q + Li hypernucleus is 26.4 MeV and that for the & + Be hypernucleus is 25.1 MeV, the cores of which (i.e., 9 He and 
w Li respectively) are known to be unbound. BWMH predicts the 6 + + 12 C binding energy to be 23.24 MeV where 
as, the same for the A + 12 C binding energy is 12.65 MeV. This shows that compared to the A hyperons, the 8 + 
hyperons will be more strongly bound in a nucleus. This finding is in consonance with the prediction of Zhong et al. 
[18] and QMF calculation [19]. 

It is pertinent to note that unlike the previous two mass formulae [3], [20] the BWMH mass formula does not 
jeopardize the normal nuclear matter properties. Fig. 4 shows that the BWMH predicts the line of stability quite well 
while the mass formula of Dover et al. [3] shows mismatch for medium and heavy nuclei. The mass formula of Levai 
et al. [20] although reproduces the line of stability, it does not reproduce the binding energy versus mass number (A) 
plot (Fig. 5a) of non-strange light nuclei such as, e Li [37]. Mass formula of Dover et al. also shows marked deviations 
for isotopes with larger neutron numbers. The sharp oscillations in the experimental neutron separation energy (S n ) 
versus A plots (Fig. 5b) are reproduced by BWMH but, not by the other two as the pairing term is altogether absent 
in the mass formula of Dover et al. [3] and, in case of the mass formula of Levai et al. [20] it is too small. The 
incompressibility of infinite nuclear matter [38] obtained using the energy co-efficient a v of the volume term is about 
300 MeV [39] for BWM (and BWMH) which is within the limits of experimental values while the mass formula of 
Levai et al. [20] predicts values in the range of 400-480 MeV which are too high to be realistic. The presence of 
Majorana term in the mass formula of Levai et al. (eqns. (6) and (8) of [20]) also poses a serious problem that the 
binding energy per nucleon diverges as mass number A goes to infinity. This violates the nuclear saturation properties. 
As BWMH is not plagued with such divergence problems and nuclear saturation properties are well preserved for 
large A, this mass formula will be useful for extension to astrophysics related problems, like equation of state of a 
neutron star. 

In summary, a single mass formula (BWMH) valid for both non-strange normal nuclei and strange hypernuclci is 
prescribed by introducing hyperon mass and strangeness dependent SU(6) symmetry breaking terms in the modified- 
Bethe-Weizsacker mass formula (BWM). The BWMH preserves the normal nuclear matter properties. Due to delicate 
balance amongst the surface, strangeness and Coulomb contributions, a maximum in q = +1 hyperon separation 
energies is found near lower A values along with a decreasing trend at larger A. This feature is absent in hypernuclei 
with q = and q = — 1 hyperon(s). Hyperon separation energies calculated by BWMH are in good agreement with 
the available experimental data on 'bound' nuclei, including |-ffe. Heavier Sigma-hypernuclei were predicted to be 
unbound for repulsive Sigma-nucleus potential [6] . As BWMH does not account for the repulsive potential, it predicts 
high binding energies for the heavy Sigma hypernuclei which contradicts the present wisdom. On the other hand, 
several authors have predicted that Theta-hypernuclei would be bound more strongly than Lambda hypernuclei. 
Quark mean field estimates for the Theta-separation energies are now available for medium to heavy-mass Theta- 
hypernuclei [19] and those predictions are in close agreement with the BWMH predictions. It is noteworthy that 
BWMH can predict separation energies of both light and heavy hyper-nuclei. In view of the proposed search of the 
elusive + hypernuclei and for many other hypernuclei which are predicted to be 'bound' (but their binding energies 
are unknown) the present mass formula is expected to provide a useful guideline. 

We gratefully acknowledge G. Levi for sending us an useful compilation of some experimental data on hypernuclei. 
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FIG. 1. Hyperon separation energy Sy versus mass number A for (a) single A (my = 1115.63 MeV) predicted by BWMH 
(solid), Dover et al. [3] (top most, dashed-line) and Levai et al. [20] (dotted) and experimental values [1, 2] (rombus with error 
bars) and, (b) double A predicted by BWMH (solid), and Levai et al. [20] (dotted) and experimental values [1, 2] (rombus 
with error bars) . In the inset the three data points of double A are shown with predictions of BWMH. In all the figures, lines 
are added only as guides to the eyes. 
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FIG. 2. Hyperon separation energy Sy versus mass number A for (a) single H (my = 1321.32 MeV) and experimental values 
[1] and, (b) single 6> + (my = 1540 MeV) separation and quark mean field calculations of Shen and Toki [19]. 
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FIG. 5. Predictions of BWMH, Dover et al. [3] and Levai et al. [20] and the experimental data [37] for (a) binding energy 
versus mass number A for Z=3 and (b) the one neutron separation energy versus A for Z—38. 
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